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Initial-state nuclear effects in proton-nucleus collisions 
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Abstract. Two important initial-state nuclear effects in hadron-nucleus collisions are considered. The ratios 
of inclusive differential cross sections for Drell-Yan dimuon production are calculated. The calculated results 
are compared to the E866 data. It is shown that consideration of multiple soft rescatterings of incident 
quarks in nuclei and initial-state quark energy loss effects allow to get a good agreement between the 
calculated results and the experimental data. 

PACS. 13.85.Qk Inclusive production with identified leptons, photons, or other nonhadronic particles - 
24.10.Lx Monte Carlo simulations - 25.40.Ve Other reactions above meson production threshold - 25.75.-q 
Relativistic heavy-ion collisions 



1 Introduction 

Nowadays there is a great interest among the widest cir- 
cles of physicists in non-trivial effects of relativistic nu- 
clear physics such as anomalous nuclear dependence in 
processes with large transverse momentum px ■ This prob- 
lem has become especially important in connection with 
the recent data from RHIC HKKIEI 

It was first observed in 1975 |Sj that high-px hadrons 
in proton-nucleus (pA) collisions are produced copiously 
in the range of px ^ 2 GeV/c. This effect ("Cronin ef- 
fect" ) was observed in fixed target pA collisions at energies 
200, 300 and 400 GeV. Cronin effect demonstrates that 
a hadron-nucleus (hA) collision can not be presented as 
a simple superposition of hadron-nucleon (hN) collisions. 
Analogous behavior was observed in collisons of heavy nu- 
clei (Pb+Pb, Pb+Au) at v /5i7N = 17 GeV (SPS, CERN) 
J2f. But the recent experimental data from RHIC showed 
strong suppression of produced hadrons in central Au+Au 
collisions at ^/snn = 130, 200 GeV 

The anomalous A-dependence can be affected by initial- 
and final-state effects, i.e. the effects before and after hard 
scattering respectively. The investigation of the final state 
could give the information on the properties of produced 
medium. But this information can be extracted from data 
only when the initial state can be reliably predicted. A 
unique tool for studying the initial state is the Drell-Yan 
(DY) lepton-pair production 0, which provides the pos- 
sibility of probing the propagation of partons through nu- 
clear matter in its ground state, with produced lepton pair 
carrying away the desired information about the initial 
state, without being affected by the produced medium. 
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In this paper two important initial-state effects which 
take place in hA collisions are studied: multiple soft rescat- 
terings of quarks of the incident hadron in nuclei and en- 
ergy loss of fast quarks in nuclear matter. In order to simu- 
late the mentioned effects in DY lepton-pair production in 
proton-nucleus and nucleus-nucleus (AA) collisons we de- 
veloped a new Monte Carlo (MC) event generator HARD- 
PING (Hard Probe Interaction Generator). It is based on 
HIJING generator S, an extension of PYTHIA for 
hadron collisions on jet production in nuclear collisions. 

Multiple soft rescatterings and energy loss of quarks 
were taken into account according to refs. |10lll| . Ob- 
tained results were compaired with E866 Collaboration 
(FNAL, USA) data at 800 GeV £2]. 



2 Studying the initial state 

2.1 Nuclear shadowing versus energy loss 

Measurements of nuclear structure functions in deeply in- 
elastic lepton-nucleus scattering (DIS) [1811141115] indicate 
clearly that parton distributions of bound nucleons are 
different from those of free nucleons. 

It is very important to disentangle between the effects 
of shadowing and energy loss, since they are similar in 
many respects. In order to describe the modification of 
parton distributions in nucleus, a variety of approaches to 
this question exist in literature [T?)llT7j . 

The first DY data suitable for such an analysis were 
obtained in Fermilab E772/E866 experiments. An analysis 
of the E772 data ^Hj was made in ref. ^Hj) ignoring shad- 
owing. A better analysis was performed by the E866 col- 
laboration using the E866 data 12 . The E866 experiment 
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extends kinematic coverage of the previous E772 experi- 
ment, which significantly increases its sensitivity to parton 
energy loss and shadowing. Vasiliev et al. attempted to 
improve the analysis of ref. |18| by including shadowing. 
However the procedure employed by Vasiliev et al. used 
the EKS shadowing parametrization |20| . which included 
the E772 DY data, which are subject to corrections for 
energy loss. Thus the EKS "shadowing" already includes 
corrections for energy loss. And this is why the analysis 
of the E866 data performed by the E866 collaboration re- 
sulted in zero energy loss. 

In 2001, Hirai, Komano and Miyama (HKM) [2H pro- 
posed nuclear parton distributions, which were obtained 
by quadratic and cubic type analysis, and determined by a 
\ 2 global analysis of existing experimental data on nuclear 
structure functions without including the proton-nucleus 
DY process. 

In our present analysis we use the HKM nuclear shad- 
owing parametrization. Since the HKM fit did not include 
Drell-Yan data, we expect to find energy loss and shad- 
owing corrections which are unambiguously separated. 



2.2 Multiple soft rescatterings of quarks 

First of all, the mechanism of multiple interactions signif- 
icantly changes with energy. At low energies a high trans- 
verse momentum parton is produced off different nucleons 
incoherently while at high energies it becomes a coherent 
process. This is controlled by the coherence length [22j 



lr. 



(1) 



where fcx is the transverse momentum of the parton pro- 
duced at mid rapidity and then hadronizing into the de- 
tected hadron with transverse momentum pt- 

For a coherence length which is shorter than the typ- 
ical intcrnucleon separation, the projectile interacts inco- 
herently with individual nucleons. The energy range of 
the E866 experiment corresponds to the regime of short 
coherence lengths. Hence effects of coherence are not so 
important here. Thus we are not going to consider coher- 
ence effects in our present analysis. 

In hA interactions a quark of the incoming hadron can 
undergo soft collisions (with small momentum transfers: 
\t\ < 1 GeV 2 ) as well as hard ones (with large momentum 
transfers: \t\ > 1 GeV 2 ) inside the nucleus. It was first 
shown by Levin and Ryskin |2.'ij that the observed Cronin 
effect can not be explained by only hard collisions taken 
into account. One also have to consider soft rescatterings 
of additive quarks of the incident hadron before the hard 
process |2*3| . 

Such a picture was suggested in ref. ^U] in the frame- 
work of the additive quark model |25l . In this approach 
the dynamics of hA interactions could be visualized in 
the following manner. Each constituent quark of the in- 
cident hadron (valon) is scattered independently of the 
other quarks (as in the additive quark model) several times 
softly, i.e. with small momentum transfer, by the nucleons 



of the target nucleus. Then a (anti-) quark-parton, which 
belongs to this quark, undergoes a hard collision with an 
antiquark (quark) of the nucleus producing an observed 
DY pair. All these soft rescatterings affect the pt distri- 
bution of the quarks of the incident hadron, and therefore, 
the pt spectrum of observed DY pairs. 

According to ref. jTH] the probability for a quark to 
undergo n soft rescatterings in nuclear matter is: 



oo 

Pn = T^T)! j dZ j dM^-(M)] 
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x p(b,z)e- aT - {b ' z} 
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Here b is the impact parameter of hA collision, a is 
the cross section of the inelastic soft quark-nucleon inter- 
action, p(r) is the nuclear density normalized to unity 
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T_(b, z) is the profile function: 



T_(M)= / p(b,z')dz 



Transverse momentum fcx distribution for the quark 
suffered n soft rescatterings can be written in the following 
form 26 : 



/n n 
nd 2 p Tl / q (PT S )<5 2 (fc T -5>Tj , 
i=\ i=l 



(3) 



where p r r i is the transverse momentum of the quark in i- 
th rescattering process. Here f q (pr) is the probability for 
the quark q to have the transverse momentum after a 
single quark-nucleon interaction: 



1 da 



c d pt 



(4) 



Finally, eq. J5J can be presented in the following form |27] : 
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where m = n — 1, K m (y) is the McDonald function of 
m order, r(a) is the gamma function, B = 2/{ky), where 
(ky) is the mean momentum of the quark in the projectile 
proton. 

The modification of primordial pt due to soft rescat- 
terings of the quark undergoing the hard process was im- 
plemented in HIJING according to the mentioned distri- 
butions (eqs. © and ©). Fig. □ represents HARDPING 
simulation results for the ratios of the DY pairs produc- 
tion cross sections in pA collisions versus pt of produced 
pairs. The ratios were calculated for different (fcy) values 
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Fig. 1. Ratio of the cross sections for Drell-Yan events ver- 
sus pt with quarks multiple rescatterings taken into account. 
Ratios calculated for a — 10 mb and for different (fey) values. 



and for a = 10 mb. Also the simulation results, obtained 
with the standard HIJING multiple scatterings algorithm 
(HIJING "pT-kick" ) |Sj as well as E866 data are presented. 

The ratio o"w/&Be is a fraction of inclusive differen- 
tial cross sections, normalized to the corresponding atomic 
number: 



cw/cBe 



1 dqP w \ 
A w dp T J , 



1 der pBc 
^Bc APT 



(6) 



where A-w and As e are atomic numbers of corresponding 
nuclei, da pW /dpx and dc pBe /dpT are inclusive differential 
cross sections for DY pairs production in corresponding 
reactions. 

Fig- 121 represents the same ratio for different values of 
quark-nucleon cross section and for {ky) = 0.4 GeV/c. 

As seen from figs. ^ and [21 considering the multiple 
soft quark rescatterings effect allows to improve agree- 
ment with the experimental data comparing to the results 
obtained without taking into account this effect as well 
as to the HIJING "px-kick" results. Also one can see the 
results strongly depend on the values of a and (ky). Un- 
fortunately, in spite of the obvious improvement of agree- 
ment between experimental and simulated data comparing 
to the original HIJING, there is a large quantitative in- 
consistency. In fig. |3| one can see the best result obtained 
within this model, where the value of p^ as 2.2 GeV/c, 
corresponding to the maximum of the ratio, is about 1.2 
times less than the experimental one. Moreover, the value 
of a = 5 mb corresponding to this figure is about 2 times 
less than the one predicted by the additive quark model 
(a w | ctnn ~ 10 mb, where cnn ~ 30 mb). So we need 
more accurate model. 

The above treatment of multiple scatterings of quarks 
was oversimplified because there was no distinction be- 




Fig. 2. The same as in fig. Q but for different a and (fey) 
0.4 GeV/c. 
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Fig. 3. The same as in fig. Qbut for a = 5 mb and (fey) 
0.4 GeV/c. 



tween constituent quarks (valons) which undergo soft scat- 
terings and point-like partons which undergo hard pro- 
cess. In fact, to obtain distribution function of the parton, 
which undergoes the hard process, of the incident hadron 
one should calculate the convolution of the distribution 
function of a valon after n multiple soft collisions G™(fcx) 
and the distribution function of the parton inside the valon 
F v (k T ): 



-Pn(^t) = G™ <X> F\ 



r 

/ d 2 p Tl d 2 pT 2 Gq(p Tl )F v (pT 2 ) 

x S 2 (k T -p Tl -pt 2 ) , (7) 
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The distribution function of the parton inside the valon 
was taken in the following form: 



Fy(k T ) = ^-e 

Z7T 



Bki 



(8) 



where B — 2/(fc q lard ), and (fcq ard ) is the mean momentum 
of the parton in the valon. 

The HARDPING calculation results with the above 
treatment of multiple interactions are presented in figs. 01 
and for different values of (fcq ard ) and (ky). All calcu- 
lations were performed for a — 10 mb, which corresponds 
to the additive quark model. One can see that the ra- 
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Fig. 4. Ratio of the cross sections for Drell-Yan events versus 
Pt calculated for different (fcq ard ) and (ky) = 0.4 GeV/c (see 
the text). 



tio (Jw/cbo is rather sensible to the variations of (fcq ard ) 
and changes weakly with the variations of (ky). The new 
model gives us much better agreement with the experi- 
ment than the previous one. One should estimate the op- 
timal model parameters as 

a = (8 - 10) mb, 
(fc^ ard ) = (0.8-1.5) GeV/c, 
(ky) = (0.2 - 0.4) GeV/c. 
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Fig. 5. The same as in fig. 21 but for different (ky) and 
(fcq ard ) = 1.3 GeV/c (see the text). 



Much of the problem originates from the impossibility 
of direct measurements of this energy loss. And the lack of 
a common agreement of the processes and mechanisms to 
be included in specifying the energy loss leads to large dif- 
ferences between the results of each separate experiment. 
The energy loss in the initial state are usually measured 
from the data on A-dependence of the Drell-Yan pair pro- 
duction in proton- nucleus collisions, see e.g. ref. ^2]- I n 
this paper we adopted the model for initial-state quark 
energy loss described in ref. [TT] . 

It is usually assumed that quark propagates from the 
surface of the nucleus to the point where the DY pair is 
produced, which would mean that the mean quark path in 
the nucleus would be (L) ss ZR A /A. But, as shown in ref. 
|29| . this value should be shortened by at least the mean 
free path of a proton in a nucleus, ~ 2 fen. This would 
substantially reduce (L) by a factor of two or more, so that 
the mean path between the point of DY pair production 
and the first inelastic interaction is actually shorter than 
the maximum possible distance to the edge of the nucleus. 
Additionally, there is some probability (dominant for light 
and medium-heavy nuclei) that the incident hadron has no 
interactions prior to the point of DY pair production. In 
accordance with the above considerations, the mean quark 
path in the nucleus can be written as |llj : 



2.3 Energy loss of fast quarks in nuclear matter 

There is another important process which affects proceed- 
ing of hadron-nucleus collisions - the enery loss of fast 
quarks while travelling through nuclear matter. This pro- 
cess also takes a great interest last years, there are a lot 
of models which are trying to describe such an effect, see 
e.g. ref. but unfortunately there is no general agree- 
ment on the value of the partonic energy loss rate dE/dz 
neither in cold nor in hot nuclear matter. 



hN f r 

L) = (1 - Wb) -g- J d 2 b J dz 2PA (b,z 2 ) 



x / dzi p A (b, z\) [z 2 - zi) 



< N / dz PA (b,z) 



xexp 



(9) 
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The exponential factor requires that there is no inelastic 
interaction of the beam hadron prior to the point z\ . Wq 
is the probability of no inelastic interaction of the beam 
hadron in the nucleus prior to the DY reaction, which can 
be written as 



Wo = 



1 



A 07 



hN 



d 2 b 



1 - e" 



'T(b) 



A<7 



hN 



(10) 



The corresponding probability distribution in L is given 
by expression [TT|: 



W(L) = W S(L) + Wi(L) 



(11) 



where 



rate —dE/dz — 3 GeV/fm are presented in fig. El which is 
the DY cross section ratios for W to Be as functions of X\ 
(the light-cone momentum fraction of the incident proton 
carried by the produced DY pair) for various intervals of 
the invariant mass M of produced DY pairs. We calcu- 
lated shadowing for the mean value of mass calculated for 
each interval as -J (M 2 ). One may expect a substantial 
scale dependence of nuclear structure functions. However 
at considering ranges of x\ > 0.3 and 4 < M < 8 GeV/c 2 
and at our statistical accuracy this dependence is not seen. 
From comparison with the experimental data, it is found 
that our calculated results are in good agreement with the 
Fermilab E866 data. 



Wi(L) = Z±- J d 2 b J dz 2PA (b,z 2 ) J dz lPA (b, Zl ) 



xS(z2 — zi — L) exp 



hN 



dz px(b,z) 



(12) 



The HARDPING procedure to simulate quark energy 
loss uses two parameters: energy loss rate — dE/dz (free 
parameter), and length of a quark path L in a nucleus, 
which is calculated according to the distribution eq. (Tfl"1) . 




Fig. 6. Ratio of the cross sections for Drell-Yan events versus 
x\ calculated for optimal value of the parton energy loss ratio 
—dE/dz = 3 GeV/fm. M is the invariant mass of produced 
DY pairs. 



The results of HARDPING calculations with energy loss 



3 Conclusion 

In this work the analysis of influence of the multiple soft 
rescatterings and the energy loss of fast quarks effects on 
DY pairs production in proton-nucleus collisions were pre- 
sented. In order to perform such an analysis we developed 
a Monte Carlo event generator HARDPING, which ex- 
tends HIJING program to the mentioned effects. 

The analysis has shown that the proper treatment of 
multiple soft rescatterings of a quark of the incoming had- 
ron allows to significantly improve the agreement between 
the simulated and experimental data comparing to the 
original HIJING. The shape of the ratio ow/cBe is very 
sensitive to variations of the mean momentum of the par- 
ton inside the valon (fcq ard ) and the cross section of the 
inelastic quark-nucleon interactions er. 

The simulation of the xi-dependence of the ratio <r\v / CBe 
showed that consideration of the quark energy loss effect 
improves agreement with the experimental data. Original 
HIJING does not take into account this effect, which leads 
to contardiction with the data. 

However one should be cautious applying the results 
of this analysis to higher energies of RHIC and LHC, be- 
cause at RHIC and LHC energies coherence effects are im- 
portant. The reason is that coherence length l c becomes 
large at high energies comparing to the typical internu- 
cleon separation, hence the projectile interacts coherently 
with individual nucleons. Moreover, the initial-state en- 
ergy loss effect is not really important at the energies of 
RHIC and LHC. Thus one can disregard energy loss and 
test models of shadowing by direct comparison to data. 
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